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Abstract. We investigate how the QCD equation of state can be reconstructed 
by a continous mass distribution of non-interacting ideal components. We find 
that adjusting the mass scale as a function of the temperature leads to results 
which are conform to the quasiparticle model, but a temperature independent 
distribution also may fit lattice simulation results. We interpret this as a support 
for the quark coalescence approach to quark matter hadronization. 
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1. Distributed mass quark matter 

The quark gluon plasma, produced in the Big Bang or in high energy relativistic 
heavy ion collisions, hadronizes. We may detect and study this particular matter by 
observing signs of its collective behavior, in local equilibrium obtaining its equation 
of state. Several statistical and hydrodynamical models have been considered in the 
soft QCD sector to describe hadron spectra. 

We have developed a massive quark matter coalescence picture^] which 
determines hadron yields and transverse spectra according to branching ratios between 
concurrent hadronization channels. Partonic level models of heavy ion reactions 
also utilized the quark coalescence picture recently |31 El IS]- The seeming entropy 
reduction problem by coalescence with an associated reduction (confinement) of color 
degrees of freedom can be resolved by assuming sufficiently massive partons around 
the hadronization temperature in the precursor matter. The necessary mass scale for 
quarks is about 300 — 350 MeV and even higher (about 700 MeV) for gluons. In [2] 
the explicit gluonic degrees of freedom were neglected. 

In order to compose low hadron masses from massive constituent quarks we have 
introduced distributed mass partons into our hadronization model [6]. The low mass 
can be obtained from the convolution of distributed masses. 

This approach led us to the investigation of an ideal but distributed mass parton 
gas. In particular we study i) the equation of state (EoS) for a continuous mixture, 
ii) the consistency of the quasi-particle picture, iii) the fit of different mass spectra 
to lattice eos data, and iv) we collect arguments in favor of a mass gap. We make a 
few interesting comparisons to non-ideal gas effects: a fixed mass M = gT, scaling 
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linearly with the temperature in the quasi-particle picture, is known to lead to a 
reduced pressure compared to the ideal gas (Stefan-Boltzmann) case even at infinite 
temperature p/psb < 1 [311111]' This approach has been motivated by partially 
resummcd high-temperature pQCD calculations. A pressure reduction even at high 
temperatures also occurs in non-perturbative lattice QCD calculations. 

Slightly below T c the pressure does not vanish exactly. While earlier it was 
attributed to the finitcness of the modeled lattice, modern scaling techniques along 
the line of constant physics made us to believe that the nonzero pressure here is a real 
effect. In fact massive resonance gas eos fits quite nicely this emerging part in the 
temperature range (0.8 — 1.1)T C [Redlich,Tawflik]. We re-evaluate these phenomena 
in the light of the distributed parton mass model. 

The particle spectrum is given as a convolution integral of the spectral function 
p(s) and a statistical factor /. The former may be characterized by dynamical factors, 
as e.g. a mass scale M, the latter by the temperature T and chemical potential 
p which are characteristic to the medium. A quasiparticle is described by a delta 
function with an arbitrary dispersion relation E p : p(E,p) = ^-5(E—E p ). The 
distributed mass parton is equivalent to a continuous, finite-width spectral function: 
p(s) = fdmw(m)5(m 2 — s) = ■ An ansatz for w(m) is equivalent to an ansatz for 

p(s). We note that field theoretical in- medium spectral functions break the Lorentz 
covariance and show a separate dependence on the energy E and momentum p. 



2. Consistent equation of state with mass distribution 

Thcrmodynamical consistency of the quasiparticle picture imposes further constraints 
on the mass distribution, w(m). This can best be seen when starting with a 
homogeneous equation of state, given as a continuous sum of partial pressure 
contributions supported by a mean field part: 

P(T,V)= / dmw(m)p m (T,p,) - $(T,/x). (1) 



Here the partial contributions are given by the ideal gas formula at a fixed mass, 

-dE k fc 3 (E k -p 
2tt 2 3 V T 



p m {T,p) = /^|^/(^_^) (2) 



with f(x) either the Fermi or the Bose (or approximately the Boltzmann-Gibbs) 
distribution and with a particular, m-dependent dispersion relation for a free particle 
fc = \/ E 2 — m 2 . The entropy density and number density, as respective derivatives 
of the pressure now contain extra terms due to the T- and ^-dependence of the mass 
distribution and of <&, resulting in the following energy density, e = Ts + pn — p: 

f , , n / n f , dw dw ( 9$ <9$\ 

e= Jdmw(m) e m (T, fi) + J dm T— p m + p— p m + I <P - T— - p,— \ . (3) 

The quasiparticle consistency requires, that the total energy is also a sum of the 
respective individual contributions plus the mean field contribution, <£>. Therefore the 
mass distribution w(m) has to satisfy nontrivial constraints 

, dw (9$ f dw d$ 

dm—p m {T,p) = — , J dm — Pm (T,p) = — . (4) 
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In order to obtain $(T, fx) the integrability condition, dj. = <9^ T <I ) , has to be 
satisfied. This, assuming that w(m; T, fi) is integrable, leads to 

dm^n m (T,fi) = J dm^- s m (T, fi). (5) 

The trivial way to satisfy this is to use a T- and /^-independent mass distribution, 
w(m). In this case the mean field part, <&, is also T and /^-independent and reduces 
to an old-fashioned bag constant. The next step is the reconstruction of the equation 
of state. The pressure (and energy density) modification ^(T, fi) is obtained from 
integrating the constraint equations J3J|. The total pressure (eqQ becomes 

p(T,fj,) = J dT J w(m,T,fi)s m (T,fS)dm— $(T c ,/z). (6) 

The mean field term, $, cancels in the combination of e + p. This fact can help one 
to guess an appropriate mass distribution consistent with lattice eos data, as well as 
with a quasiparticle picture. 

We consider from now on a particular class of mass distributions, which depend 
on the thermodynamical environment parameters T and \x only through a single mass 
scale, M(T,fi): 

Tfl 

w ( m ,T,p) = — /( — ). (7) 

The normalization integral for w is inherited by the shape (form factor) function f(t): 

J w m dm = J f(t)dt = 1. (8) 

All medium dependencies are concentrated on M(T, /<), which should satisfy the 
constraint stemming from the integrability condition eq.(0. 

Now we try to guess the proper mass distribution in order to arrive at a pressure 
resembling lattice QCD simulation results. A particular, analytically integrable ansatz 
for the mass distribution shape is given by 

f(t) = Aexp(-at-b/t). (9) 

Its normalization can be obtained from the general formula 

dt t P- 1 e - at - b / t = 2(b/a) p / 2 K p (2Va~b). (10) 

The pressure contribution normalized by the Stefan-Boltzmann pressure, psb — Po, 

a(g) = f w(m)^dm - [°° dt f(t)^K 2 (gt), (11) 

J PSB Jo 2 

with g = M/T, is given by 

^3 r , 

<j(g)=g 2 ^ [K2(Vb(Va^+Va~^))K 2 (Vb(Va~T^-Va~g))\ ■ (12) 

In figure ^ we show the entropy density normalized by the Stephan-Boltzmann 
value (belonging to the massless ideal gas), s/ssb = (e +p)/(e + p)sb as function of 
the temperature (a) and the corresponding mass distributions (b). 
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T / M(T,0) m / M(T,0) 

Figure 1. (Color online) Normalized entropy - temperature curves (left) for 
different parameters of the exp(— am/T c — bT c /m) mass distributions (right) 
compared with lattice QCD results^ (boxes left). 

3. Fits to lattice QCD results 

The lattice data of p/psb can be related to a(g) by fitting the results. Since due to 
ca. i|ll|) er(0) = f(t)dt, the normalization of the mass distribution requires <r(0) = 1. 
In this casep = Psb would be satisfied at infinite temperature (g = T c /T). As a matter 
of fact lattice data reach only p/psb = 0.8 at T w 4T C instead. The quasiparticle 
model with a fixed mass which is growing proportional to the temperature actually 
predicts such a deviation even in the infinite temperature limit: in this case cr(0) < 1. 
Therefore it is important to jugdc without prejudice whether lattice QCD data suggest 
a <t(0) value lower than one or not. 

In our understanding a plot in terms of g = T c /T mediates a better picture of the 
extrapolation to the infinite temperature point at g = 0. The lattice QCD eos data 
do not contradict to <r(0) = 1. It is of course still imaginable that a bending down 
occurs at high values of T (low values of g) not simulated so far. The tendency with 
growing lattice size (i.e. the comparison of N t = 4 and N t = 6 data), however, seems 
to support er(0) = 1, a tautologic consequence of the distributed mass model. While 
this high-temperature behavior is well fitted by the pure exponential a(g) function, 
the part below T c cannot be recovered this way. A more steeply rising trial function 
is needed. 

Figure [3 plots the expectation value of the mass as a function of the temperature 
for the mass distribution © obtained by adjusting the pressure to the lattice result. 

Figure |2| presents the fitted pressure to lattice data and the interaction measure, 
(e — 3p)/psB- The mass distribution Q as a function of the mass m and the 
temperature T is plotted in figure 0] 

At the end of this section we consider the mass distribution © and investigate 
the temperature and baryochemical potential dependence of its mean mass scale 
parameter, M(T,fx). 

The explicit showing of the other parameter(s), i.e. &, we shall suppress in the 
followings: by plotting the expectation value of the mass instead of M this is well 
founded. The entropy density becomes in this special case 

dp A 3 / n f M ,,dM\ A cM> 
S= df = V* T {^-T^+y-T^Wj^J-W (13) 
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Figure 2. (Color online) The expectation value of the constituent mass as a 
function of the temperature as it follows from adjusting the mass distribution 
width parameter to lattice equation of state data. 
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Figure 3. (Color online) The normalized pressure (left) and the interaction 
measure (e - 3p) /T 4 (right) obtained by using the mass distribution with 
adjusted mass sacle M(T). 



The quasiparticle consistency is equivalent to the principle that all terms, which were 
not there for a constant-mass, no-mean-field calculation, cancel: 

<9$ 



tt 2 dT 



7<7 = 



dT' 



(14) 



leaving us with 



s = ~2 ( 4cr T ~ f?°V ~ ■9T cr ') ( 15 ) 

Similarly the number density (here a conserved number density, in QGP the one third 
of the baryon charge density, i.e. the quark minus antiquark density) becomes 



d P A r r3 I 



dM 



<9$ 

The second quasiparticle consistency equation reads as 

9$ 



-A T3 dM 

7T 2 (9/1 



7<7 



d/j, ' 



(16) 



(17) 
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Figure 4. (Color online) The mass distribution for different temperatures, 
w(m,T) = f(m/M{T))/M(T), with an adjuted mass scale to lattice QCD 
pressure results. It has allover a finite width, its local maximum follows the 
behavior of temperature dependent mass known from the quasiparticle model. 
T c = 165 MeV has been used. 



leaving us with 



(18) 



The energy density is combined to be 

e = Ts + un-p = -^T 4 (3a - ga'\ 7 + (19) 

It is noteworthy (although known for long) that the combination e + p is independent 
of the mean field part: 

e + v A , , 
"j/ = ^7 (4<r -<?</) (20) 

Now we investigate the integrability condition in this special case. The /i-derivative of 
the T-consistency equation <|14[1 (we shall call left hand side, LHS) has to be equal to 
the T-derivative of the /i-consistency eq. (|17|l . Omitting two common terms in both, 
namely ja'd 2 M/dTdn and fff^Tff", we are left with 



dT dfi 

_ A_ 2 , ,dM_ 
dfldT ~ ^ 1<J ~dT ' 



9 9$ A 2 ( ^ M „\ dM 



from the equality of the above expressions it follows 

, .dM / ., u 1 , „ i\ dM n 
° i^F + 2 V + £yV - 3 7 a' -^- = 0. (22) 



8T ' V" T 1 J dfx 

In this partial differential equation for M(fx,T) a depends only on g = M/T and 7 
depends only on a = \ijT. The solution can be characterized by M(T, ^i)=constant 
lines on the T — fi plane. Such curves follow an ordinary differential equation obtained 
from (22) : 

du 7 / cr" \ u 
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The solution has the form [i = Ta{g) whence we obtain 



a 



9-7 - 3- (24) 



1 g?7 / dfi /i 
7 da \ dT T 
The expression in the bracket above reduces to 
dfi, fj, da da 

dT-T =T dT=-% (25) 
while gT — M is constant. Therefore we obtain the following differential equation: 
g dj da a" 

-~-T-r=9 3, 26 

7 da ag a 

which becomes an integrable problem for 7(3): 
d n a" 3 

_ — ln 7 = — --. 27 
dg a 1 g 

Its solution is given by 

7 = Kg 3 /<r'(g). (28) 

This, replacing the definition of 7 in terms of fi/T, and g = AI/T with constant 
M = mo one obtains a definite integral curve starting at T = Tq for /1 = 0. That 
means that the constant K is related to To. By inverting the cosh function the 
following explicit solution emerges for the M — mo constant lines: 



fi = T In 

with 



(z{T) + y/ z *(T) l) (29) 



z{T ) =11 + ^ T i a '^ 0/T ^ (30) 



2 7Q y r 3 a'(mo/T) 2 7Q ' 

The temperature dependence of the mass scale at fi = 0, M(T, 0) = mo, has to be 
inverted in order to obtain Xo( m o)- This positions the crossings of M = mo constant 
characteristics with the temperature axis. We note that in every model with baryon 
- antibaryon symmetry these lines have zero ^-derivatives at (i = 0. 



4. Arguments for a mass gap 

The o{g) function in the Boltzmann approximation is given by This is a so called 
Meijer K-transform (a generalized Laplace transform) which can be inverted by 

/(*) = ~ r^^^dg. (31) 

This presents a peculiar problem, whether there exists a unique mass distribution fit), 
with the mass scale parameter kept temperature and chemical potential independent, 
to any a(g) function extracted from an equation of state (e.g. from lattice QCD 
calculations). The shape of such a mass distribution is not arbitrary. We shall explore 
this possibility in a future work. 

Utilizing the parametric integral one easily derives a useful relation between 
moments of the mass distribution, f(t) and the eos fit. We obtain 

1 T (s±l) 

M n = (t-' n ) f = -= . V 2 . ' - (g"' 1 )^ = cnig 71 - 1 ),. (32) 
\ I j 0fr f + 2 Tn W ' W ' K ' 
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Figure 5. Continuation of the medium dependence of the single mass scale to 
finite chemical potential, M(T, fi), using a 12:16 light quark - gluon mixture with 
mass distribution. It have been obtained by solving the integrability condition 
1221 starting from the temperature dependent scale M(T) fitted to lattice QCD 
results at zero chemical potential. The mo = M(T, 0) values are given as 
m = 0.341,0.129,0.146,0.184,0.226 and 0.270 GeV from above. 



A roughly approximate, qualitatively correct fit to the lattice eos data is represented 
by the straight line, <j(g) = 1 — g/g c , with g c ss 1.15. All the n > moments of this 
expression arc finite, meaning that the inverse mass moments are also finite: 



This is possible only if the mass distribution has a finite mass gap, or it approaches 
zero more than any polynomial in the inverse mass 1/t. Both possibilities represent 
an interesting spin off of the eos studies. The sizeable reduction of the pressure at low 
temperature (large g), which causes the n > moments of a(g) be finite in general, 
is related to confinement. This requires finite moments of the inverse mass with the 
mass distribution, which is possible only with a corresponding suppression of the low 
mass part, alike exp(— M/m), or with a mass gap. 
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